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INTRODUCTION 


Nature of the Problem 

It is well-known that aircraft can be greatly affected by cross- 
winds and gusts during take-off or landing. 

The situation is particularly critical for vertical and/or short 
take-off and landing (V/STOL) aircraft which may be operating in con- 
gested urban locations and thus subjected to the wakes of hangars or 
tall buildings. 

Wakes are regions in which the mean velocity is less but the 
turbulence is greater than that in the unaffected flow. Semi steady 
state vortices may be a feature of tie wake in the near vicinity of 
the structure. 

Aside from these general statements, little is known about the 
wakes of real structures in natural atmospheric flow. It follows that 
potentially hazardous situations for aircraft operations conceivably 
could be avoided if there was greater knowledge of the structure of 
wakes and how that structure is related to obstacle size and shape, and 
ambient wind. 

Past Studies 

Frost (1973) has provided a very useful review of over 50 papers 
and symposia proceedings dealing with air flow over bluff obstacles. 

His review reveals that studies of models in low- turbulence wind 
tunnels constitute by far the greater portion of the literature on the 
subject. These studies deal with flow over plates, forward or rear 
facing steps, or buildings. 
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Typical of the sort of flow found over buildings or plates is that 
given by Halitsky (1968) and shown here as Figure 1. 



Figure 1. - Flow zones near a sharp-edged building (Halitsky, 1968). 


It is seen that flow separation occurs near the forward edge of the 
building roof and that reattachment occurs at some distance downstream 
within the wake. A closed circulation is thus delineated which has 
been variously called, a cavity, a bubble, or an eddy. 

Evans (1957) has concentrated on determining the dimensions of 
this feature of the wake in a wind tunnel study of buildings of 
various shapes and roof types. In addition, he shows that the flow 
pattern in the horizontal plane not only depends upon the shape of 
the structure but varies considerably with the orientation with 
respect to the wind. Unfortunately, he made no attempt to develop 
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analytical relationships in this Study. 

Hal its ky provides equations dealing with the wake dimensions and 
mean longitudinal velocity distribution, but they are applicable to a 
suspended plate and therefore of limited value for the natural setting. 

A very elegant theoretical trea^m^nt for the flow in the vertical plane 
is given by Couni han et al. (1974). 

'me numerical simulation studies of flow around buildings (Hirt 
and Cook, 1972; Hotchkiss, 1972; Thomas, 1971; Frost, et al., 1973) 
show results similar to those of the wind tunnel studies. 

It is generally agreed that v/akes in the natural setting may 
differ considerably from those obtained in the wind tunnel or by 
numerical simulation. This is because the real wind in passing over 
the surface boundary is itself extremely turbulent and subject to random 
variations in speed and direction. Unfortunately, little v;ork has been 
done in the natural setting. Sexton (1971) describes measurements 
made on the upwind side of a building which show a downwash vortex 
(see Figure 1) extending about one building height upstream. Measure- 
ments of average wind speed at the 3-m level above ground show 
reasonable agreement with wind tunnel data. 

A measurement program in the wake of an airport hangar has been 
reported by Colmer (1971). Instrumented masts were placed along a 
line perpendicular to and through the middle of the long side of the 
hangar at distances of 5, 14, and 23 building heights (10 m) . Another 
instrumented tower 190 m upstream provided reference data. Measure- 
ments made during a single 40-min run give indications of the extent 
of the wake, .-the intensity of the turbulence, and gradients of the 
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mean flow within the wake. 

Purposes of Stud y 

The present study was intended as a feasibility experiment to 
determine the extent to which a research aircraft could be used as a 
rapid- turnaround, data-acquisition system to yield useful data for 
aeronautical safety applications. Such a system could be used to probe 
the wake from many buildings and obstacles and would avoid the labori- 
ous and time-consuming method of measurements obtained from towers. 

It is evident that the starting point in such an experiment would 
be the establishment of standard data against which to compare the air- 
craft data, li .'s can be done only by making measurements in a wake 
with instrumented towers and developing a model based upon the perti- 
nent parameters of the flow. The next problem would be the determina- 
tion of some readily-obtainable parameters in the response of the air- 
craft which indicate the likely structure of the v;ake. 

Approach to Study 

One obvious signature of wake flow v;ould be the greatly enhanced 
turbulence at its boundaries and interior. Studies have been published 
of the response of an aircraft to a field of turbulence, e.g., Jones 
(1969) and Kaynes (1971). However, for this to be a fruitful approach, 
it is required that the turbulence be homogeneous and steady state. 

Furthermore, as Colmer (op. cit.) points out, a conventional 
aircraft at take-off or landing speed would pass through the wake of 
an isolated structure in one to two sec. The response of the aircraft 
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would not be to the energy of the turbulence at various frequencies, 
but rather to the gradient of mean wind speed across the wake. A VTOL 
aircraft could be affected by the turbulence but the gradient of mean 
windspeed in the wake again is expected to be the dominant factor. 

Therefore, the main thrust of this study has been to obtain mean 
values of the flow. The intention was to use these measurements in 
conjunction with theoretical relationships developed by means of 
dimensional analysis to establish a model of the flow in the wake. 

The flight program v/as never accomplished. This was due in part 
because of down-time of the aircraft, in part because of delays in 
getting the aircraft instrumented, but mainly because the appropriate 
meteorological conditions (east winds) did not occur at times when 
the aircraft was operational and available. 
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EXPERIMENT DETAILS 


Study Site 

The experimental site was a large hangar (100 ft x 123 ft x <10 ft) 
at what had been a U.S.A.F. airbase at Bryan, Texas, but which now is 
a property of Texas A&M University. 

Figure 2 shows the site. The area to west of the hangar is a 
concrete apron to a distance of about 600 ft. The area surrounding the 
hangar on the north, east and south sides is grassy field. Some one- 
story buildings lie 330 ft to the north of the hangar and two more are 
shown to the south-southeast. Also shown on the north side of the 
building are a small shed, a trailer truck and what had been the flight 
control tower, a steel girder structure about 70 ft high. 

The climatology of the area is such that the prevailing wind is 
southerly from late spring through early fall. However, the most open 
and suitable area for setting up instrument^ and carrying out flight 
measuremenus is the apron. Therefore, an array of spots in a polar co- 
ordinate grid were painted on the concrete. These were numbered and 
lettered as shown in Figure 2. Favorable wind directions for experi- 
mentation would be those from northeast through southeast, conditions 
which prevail infrequently, but occasionally after a cold-front passage. 

A second grid was laid out on the south side of the building, as 
shown in Figure 2, to permit measurements during northerly flow. This 
condition is more common after cold-front passages in the area. Flight 
measurements in this area were ruled out because of power lines which 
would intersect any desirable flight path. 




Figure 2 - Experimental site. Points represent grids used to locate 
instruments. Points are 40 ft apart in south grid. 
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Ground Measur ements 

The instrumentation employed was the Climet system for measurement 
of wind speed (model 011-1) and direction (model 012-1). Gill propeller 
anemometers (model 27100) were used to measure vertical motion. The 
instruments were mounted either on portable tripods for measurements 
up to 10 ft above the ground or on three large towers each with a heavy 
base mounted on wheels. These towers could be moved with a small 
tractor over even ground; furthermore, winches on these towers per- 
mitted instrument placement at any height from about 5 ft to 60 ft 
above ground. The unevenness and softness of the ground made it im- 
possible to use these towers south of the hangar. 

An earlier project was carried out in 1973-74. During that 
phase of the work one of the towers v/as permanently located at a 
point 150 ft to the east-northeast of the hangar (see Tower 1 in 
Figure 2). During 1975-76 the towers were unavailable most of the time. 
Consequently, a tripod generally was used for the upwind measurement 
and the placement varied from one experiment to another. 

All experiments were composed of a series of approximately 15-min 
runs with data measured simultaneously at four or five locations. One 
set of instruments was always upstream from the building; these 
measurements were considered to be the reference or free-stream values. 
The intention was to operate within or at the edges of the wake with 
the other towers or tripods. After each run the towers and/or tripods 
could be relocated and on the tall towers the instruments could be 
moved to a new height. Thus, with fairly steady wind for two- three 
hours, the wake could be probed in a number of places. 
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The recordinq system, analog tape units, was located in the 
northeast corner the hangar. Thus, 500-ft cable lengths permitted 
instrument placement out to about 400 ft to the west or south of the 
building. 

A number of problems were encountered with the recording of data 
on tapes and the later conversion of the data from analog to digital 
form. Consequently, many data were lost, e.g., all vertical motion 
data for the experiments run in 1975-76. It was finally decided in 
the spring of 1976 to abandon the use of the tape recorder and a 
switch v/as made to chart recorders. 
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RESULTS OF THE STUDY 


E xperimental Data 

Those data that were recorded on magnetic tape were processed 
as follows. Analog values were taken from the tapes at each tenth of 
a second and averaged to give a 1-sec value. The 1-sec values were 
transformed to digital form and stored on a second tape for further 
analyses by computer. 

The results from three experiments done in 1974 are presented in 
Tah> 1. The columns in the table show the following, going from left 
CO right: (1) the date of the experiment and the run number, (2) the 

tower numbers, (3) the location of the given tower within the grid (see 
Figure 2), (4) the height above ground of the instruments on the 
given tower, (5) -(6) the average wind speed and direction in the 
meteorological system, and (7) the average vertical motion. The 
vertical motion for Tower 2 is missing for all three experiments. A 
W •'r an S has been used as a prefix in column 3 of Table 1 to indicate 
whether the observation site was in the "west grid" or the "south grid." 

The data obtained during 1975-1976 are shown in Table 2. The 
experiment on 6-6-76 was the only one for which the west grid was 
used. For the experiments on 11-12-75 and 6-6-76, the recording of 
data was accomplished by Esterline-Angus chart recorders, as explained 
previously. 



Table 1. - 1974 Experiments 


Date 

and 

Run 

Tower 

Location 

Instr. 
Ht. (ft) 

Wi nd 

Speed (m/sec) Dir (deg) 

Vertical 
Motion (cm/sec) 

4-17-74 

1 

Fixed 

8 

3.0 

11.5 

3.7 

8 

2 

S-A3 

8 

2.4 

16.0 



3 

S-B3 

8 

2.2 

6.2 

-0.1 


4 

S-C3 

8 

2.2 

30.2 

0.8 


5 

S-D3 

8 

3.1 

22.2 

-2.0 

4-23-74 

1 

Fixed 

8 

2.5 

29.9 

-3.2 

1 

2 

W-J4 

8 

1.8 

47.0 




3 

S-A2 

8 

1 .3 

355.0 

-10.4 


4 

S-B2 

8 

1.0 

11.5 

-4.6 


5 

S-C2 

8 

1.7 

58.0 

18.3 

4-24-74 

1 

Fixed 

8 

2.0 

75.9 

-0.6 

1 

2 

W-F5 

30 

1.8 

88.8 

— 


3 

W-F5 

60 

3.2 

188.6 

-2.8 


4 

W-E5 

8 

1.7 

86.4 

-9.4 


5 

W-65 

8 

1.5 

94.6 

-1.1 



Table 2. - 1975-76 Experiments 


Date 


and 



Instr, 

Wind 

Run 

Tower 

Location 

Ht. (ft) 

Speed (m/sec) Dir (deg) 


10-16-75 

1 

1 

2 

4 

Fixed 

B1 

Cl 

CO CO CO 

6.8 

3.1 

3.1 

333.4 

351.8 

101.7 

2 

1 

Fixed 

8 

5.1 

337.3 


2 

C2 

8 

1.3 

68.0 


4 

D2 

8 

2.8 

49.5 


3 

1 

Fixed 

8 

3.6 

343.7 


2 

E3.5 

8 

3.2 

4.6 


4 

F3.5 

8 

4.1 

351.7 


5 

G3.5 

8 

4.4 

342.5 


10-17-75 


1 

1 

2 

4 

5 

Fixed 

B1 

Cl 

D1 

8 

8 

8 

8 

3.4 

1.0 

1.1 

2.8 

335.0 

346.4 

91.7 

82.2 

2 

1 

Fixed 

8 

5.2 

334.3 


2 

E2 

8 

2.2 

12.8 


4 

F2 

8 

3.8 

331.5 



Table 2. - Continued 


Date 


and 

Run 

Tower 

Location 

Instr. 
Ht. (ft) 

Wind 

Speed (m/sec) 

Dir (deg) 

3 

1 

Fixed 

8 

5.4 

317.5 


2 

E3 

8 

2.2 

16.6 


4 

F3 

8 

4.7 

340.5 

4 

1 

Fixed 

8 

6.0 

310.0 


2 

E4 

8 

2.5 

9.5 


4 

F4 

8 

5.0 

340.9 

■ 5 

1 

Fixed 

8 

6.0 

321.9 


2 

D4 

8 

3.1 

19.7 


4 

E3 

8 

4.6 

350.0 

6 

1 

Fixed 

8 

5.3 

301.3 


2 

El 

8 

2.2 

27.1 


4 

F2 

8 

4.1 

336.1 

11-12-75 

2 

1 

Fixed 

8 

8.7 

324.6 


2 

C3 

8 

4.1 

27.0 

3 

1 

Fixed 

8 

8.7 

319.1 


2 

B6 

8 

7.1 

346.4 



Table 2. - Continued 


Date 


and 

Run 

Tower 

Location 

Instr. 
Ht. (ft) 

VJind 

Speed (m/sec) 

Dir (deg) 

6-6-76 

1 

1 

Fi xed 

8 

2.9 

49.7 


2 

W-K7 

8 

2.8 

45.8 

2 

1 

Fixed 

15 

3.2 

51.5 


2 

W-K7 

8 

2.7 

61.4 

3 

1 

Fixed 

30 

3.7 

56.1 


2 

W-K7 

30 

4.2 

70.4 

4 

1 

Fixed 

60 

3.2 

48.6 


2 

W-K7 

60 

3.6 

58.9 

5 

1 

Fixed 

8 

2.4 

65.6 


2 

W-G9 

8 

2.8 

63.8 

6 

1 

Fixed 

15 

2.9 

63.6 


2 

W-G9 

15 

3.0 

71.1 

7 

1 

Fixed 

30 

3.7 

54.9 


2 

W-G9 

30 

4.1 

60.2 



Table 2. - Continued 


Date 


and 

Run 

Tower 

Location 

Instr. 
Ht. (ft) 

Wind 

Speed (m/sec) Dir (deg) 

8 

1 

Fixed 

8 

3.5 

60.9 


2 

W-H6 

8 

2.7 

69.6 

9 

1 

Fixed 

15 

3.6 

72.5 


2 

W-H6 

15 

2.6 

73.3 

10 

1 

Fixed 

30 

3.2 

57.3 


2 

W-H6 

30 

3.0 

54.9 

11 

1 

Fixed 

60 

3.8 

63.7 


2 

W-H6 

60 

3.9 

65.9 
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Theoretical Development- 


We first establish the nature of the equations which must be used 
in the analysis by a scaling procedure, along the lines of that given 
by Tennekes and Lumley in their discussion of the turbulent wake (A 
First Course 1n Turbulence, pp. 104-109). 

Accordingly, we start with the equations o? mean motion and the 
continuity equation (the coriolis terms and molecular viscosity terms 
are neglected); 


3t 3x gy 3 z p gx 


I 

p 



+ ^ (puv) 



( 1 ) 


4+ + V-51+ w^= --1-5P - 1 

3t gx gy gz p ay p 


(puv) + 1^ (pv^) ,+ ■— (pvw) 


ax 


. ( 2 ) 


4 + + w^ = - 1 

at gx gy gz p gz 


- g - 


(puw) + ~ (pvw) + |-- (pw^) 


gz 


13) 


^ + = 0 
ax ay az 


(4) 


At this point, we assume that the flow in the wake is similar to 
that shown in Figures 3 and 4. In these drawings, two scale lengths 
are introduced. The first, ii, corresponds to that defined by Tennekes 
and Lumley as the point at which U-Uq = ^ U^. The scale length, 
in the z direction follows from the discussion of Chang's wind tunnel 
work by Plate (pp. 162-165). It marks the surface of U=% Uq'. It fol- 
lows that both i and s ,2 functions of x. 
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Then, following Tennekes and Lurnley, we can say that 


f = ""''I?' ^<VL) ■ 


( 5 ) 


where L denotes a scale length in the X direction. From the continuity 
equation, the orders of magnitude are: 


li. 


V 


W 


r + ir* ti' ° 


Assume: 




and 


f\j 


L'h 

10 L 


where ) and signify order of magnitude. 

Therefore, 


V 


-ff\ 


I»m\ 




U_5-9\ 


, W = 


lOL 


( 6 ) 


We now assume that IJ ] Oq, UoX/Uq. Then, assuming 

^ 10 10 

steady state conditions and comparing the magnitudes of the terms in 
(1) gives: 


U3U 

3x 


'\;U( 




Uf 


U, 


V :s 

3y L 

W /^S*'2 

3l" m 


\U 10 lOJli 100^1 ’ 
U.^Li /U \ U 2 Uq2 


% 

u 


I L 

0 UoU 
£2 


% 


1000£i ’ 


S 1 

a. 




lOL 10 100£i ’ 


Un2 


3U^ U' 

3X L 1000£i 


; where it was assumed u 'v -|q Uq , 
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o 1 UO 

3UV ^ ^ 1 

3y Zi 100 ’ 


im r, 

dZ ^2 100 5,2 


Since V, w u 


Therefore, the approximation to (1) is 


n Ml + + 3UW 1 3P 

^ 3 X 3y 3 Z ~ p 3 X 


Equation (2) terms become: 

,, 3V ,, f^s^'A 1 _ ^qUs«'1 
37 " ^0 h- L - 7T— 


10005,1 ’ 


,3V fV'A 1 

^ 3y '^ \ L / £i 1T~ 1000)1] 


3V f^s^2\ 1 U^^£i ^ 


\ iol/ \ l / £2 ToiT" 10000 5,1 


377 u2 1 Uo^ 

3x ^ L ^ 100 L 1000)li 


3y £1 100£i ’ 


32 £2 100£2 100^1 


(Note, £1 £ 2 ). 


Then, (2) becomes: 


3\7~ ^ 3VW _ 1 3P 

3y 3 Z " P 3y 


From (3). 


U — U 


fU^£2\ 1 UoUg£2_^ Uo^ 

ol^mj L ■ lOL^ 10,000£2 
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V 

^ 3y L 


1 

10 




1 _ ~s 

" I 


U .^^2 




Uo^ 

ToToooTi' ’ 


'dz llOL / 5-2 


lOOL^ 


Uo' 




10^i^2 


iuw -X, Ml -i__ 
3X L lO^ti 


U 


?-\ U 2 

'll — P 

100 / 1000^1 


3vw „ ^ 

9y 5-1 1005-1 


^ Mi ^ .Jlil. 

9Z 5,2 1005-2 


Thus, 3vw , 3 _ 1 3P _ ( 

w ~z p -s?- 9 • ' 

We identify P' as the pressure anomaly due to the effect of the build- 
ing and Pq as hydrostatic pressure, i.e., 

P = Po + P' . 

If we ignore synoptic effects, then we can assume v^Pq~ 0. Also, if we 
assume that the mechanical effect of the building on the flow is much 

larger than any buoyancy effects, we can use the hydrostatic approxi- 

mation. Thus, 


Ik + kMl + Mi = _ i k£i 
3x 3y 3z " p 3x ’ 


3V^ + 

3VW _ 

_ 1 

3P' 

9y 

3Z 

P 

3y 

3 VW , 

3W^ _ 

_ 1 

3P' 

8y 

3Z 

P 

3Z 


(10) 

( 11 ) 


( 12 ) 
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From Chang's work (Plate, Fig. 4.10), it may be assumed that 
P' = Uo^ 

Then letting 

P' = A (x,y,z) Uo^ 


U 




gives: 


p 

dx 

3X 

P 

1 

DP' _ 

SA 

Uq2 

p 

3y 

9y 

P 

1 

OP' _ 

Ml 

1) 2 

p 

DZ 

DZ 

P 

previous 

scale ana 

A 


1 

Uq2 

L' 

p 

100 

2,1 

A 

1! 2 
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Uo^- 


P 

100 


A 

V 

1 

V 


p 

100 

~^2 


G\ 


& 


Ur=’-' 


A_ 

X.1 p 


^2 


L^o^' 

p 


P .. 

or A 'x- Too 17 ’ 


or h 'h 


or l\ 


100 

p 

Too 


It is seen that if it is reasonable to expect that 2.1 and 2.2 are 
appropriate scale lengths for pressure as well as motion, then we must 
have L'/ 2 ,i Thus, the scale length for pressure in the X direc- 

tion would have to be much smaller than that for scaling motion. On 
the other hand, it seems more reasonable to expect the scale length 
for pressure in the x direction also to be the same as that for 
velocity. Then, L' L 'v 102i. Accepting the result from the second 
and third equations of motion that A p/1 00 then implies that 

1 8P' J ^ 

p dx 1000 Zi 


9 
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which is negligible when compared ‘ to the other terms in (7). 

When the Reynold's stress terms are expressed in terms of the 
familiar mixing length concept, the system of equations becomes: 


3U 

3X “ ^ " 8Z 


g2v , 3 

3y2' 3Z 



K 3 

^ ap" az 



Kiy 


az 

1 

= i 

aP 

p 

3y 

= i 

aP 

p 

az 


ax 


+ 


9y 


= 0 


(13) 

(14) 

(15) 

(16) 


At this point it is assumed that the flow in the wake is composed 
of a basic state plus a perturbation. The basic state is taken to be 
the undisturbed, upwind value Uq. Therefore, 


U = Uq(z) + U' (x,y,z,t) 

Furthermore, it is assumed that Uo(z) can be represented by the familiar 
log law. Then, 


And, K = k U* z. As a consequence (13) becomes 


U 


au' 

oax 


K 9_ 

^ “ap “ az 



= 0 


( 17 ) 


Long (1963) and Calder (1967) have emphasized the importance of 
using the controlling equations when they are known for purposes of 
generalizing the dimensional analysis procedure. Frequently, the 
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number of dimensional ratios can be reduced below that obtained by the 
usual dimensional analysis procedure and, therefore, a more powerful 
relationship is obtained. Such is the case here and it is for this 
reason that relationships (13) - (16) were developed. 

In the present instance we use the equations of motion plus the 
"balance equation" which is derived from the former by using the con- 
tinuity equation. 

Thus , 


_ i v2 p' 

p 


= Uo 


a2u' 


(18) 


We assume the following boundary conditions. 


U' = 0 at y = ± Q AT 

LI' = -Uq = 0 at z = Zg; U' = 0 at z = H . 
U' = 0 at X = U' = -Ug at x = 0 . 

V = at y = ± Q Ar, V = 0 at z = 0, H. 

ma A 

W = 0 at z = 0, H; W = 0 at y = + Q Ai. 


P' = 0 at y = ± Q Ai and z = H . 


3P7p 3 ^ 

3Z 3Z 



at z = 0 . 


Here, the lateral boundary of the wake is assumed to be at 
y = ± Q Ai and H is the to;; of the wake. Q is a dimensionless factor. 
It is assumed the H is 3A2 at best and that A 2 = height of obstacle. 
Following Long's procedure, let 

[X] = A [U-] = D [P7p] = G [Aj = J 

[y] = B [V] = E [ k ] = h [A 2 ] = k 

[z] = C [W] = F C Ug ] = I 



24 


When these relationships are substituted into (14) - (17), (18) 
and the boundary conditions, we obtain: 

[ X ] = A C U' ] = hAC"^ [ P'/P ] = hV^ 

C y ] = B [ V ] = h3“' [ K ] = h 

[ z ] = r [ W ] = hC"^ [ Uo ] = hAC"^ 

[^ti ] = B [ £2 ] = C 

It is assumed that U', V, W, and P'/p are all functions of 
X, y, z, K, Uo, ^ 1 , H' 2 . Consequently, it is found that: 


Ui = f IJL >L 

Uo ^ \UoZ ’ z ’ y 

= ll -f /ji_ . ^ 
do xy ^ \Uoz ' z ’ 

i_ = 1_ f. . X 

IIq xz ' l^UoZ ’ z ’ 

P' _ 1 f f K X 
pUo^ ^ ** \UoZ ’ z 





’ T 


9 



9 




9 


(19) 



( 21 ) 


( 22 ) 


In the following we concentrate on finding the functional 
relationship for U‘ for conditions that match those under which the 
observations were taken. Accordingly, it is assumed that the basic 
flow, Uo, can be represented by the log profile, i.e., 

U* 

Uo = 1 ^ In z/zo 

for which K = k U*z. Then, (17) becomes 


, 9^U' 
^ 3z^ 


0 . 


kU* 3x 


3^U' _ 9U' 
9y^ ~ 9z 


(23) 
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We identify z = Z, the height of the observations so that (23) becomes, 


aU' 32 u' ^ 7 a2u' _ Q 

kU* ax " ay^ “ az “ az^ 


(24) 


Thus, the coefficients in (24) are treated as constants. 

A multitude of expressions for U' were tested to see if they 

could serve as a solution to (24) keeping (19) in mind. One ppssi- 

bility is the relation - 

- ' - U' = Uo(z) ^ A (?) , (25) 

where n = constant, x = x/h, ? = y/2-i(x) ; h is taken to be the 

height of the building. When (25) is substituted into (24) we obtain 


d2A 

dF 


"H 


U„(Z) (Uti 
kKM 



Uo(Z) 

n A = 0 

khU*Z - J 


(26) 


We wish the solution to describe a self-preserving flow. This condition 
will be met if 


/2khU'F-\ _ 1/2 

FIjFzT/ ^ 


(27) 


Then (26) becomes 

0 + ? ^ - 2n A = 0 . (28) 

We expect that A(?) should be a maximum at ? = 0 and symmetric 
with respect to ? = 0. The values of ^i(^ for the experimental con- 
ditions in this study range from about 7 ft to 15 ft. Therefore, A(?) 
should approach zero as ? goes to plus or minus five or six. A function 
which will behave in this fashion is given by B e , where B and a 
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are constants. When this function is substituted into (28), it is 
found that it will be a solution providing a = 1/2 and n = - 1/2. The 
statement for (25) then becomes, 



The value of B can be determined by setting U' equal to at 5 = 0 

on a plot of the left-hand side of (29) versus 6. 

Another possible solution can be obtained by expressing A(^) as, 

00 

/^(5) = a. 5^ , (30) 

where the coefficients a^. are to be determined by substitution into (28). 
The symmetry condition requires that the final expression shall contain 
no odd powers of f. Accordingly, it is found that 

A(^) = ao + X (3i) 

i=l 


where ao can be determined in terms of by setting 5=0 when (31) 
is substituted into (25). The difficulty with this solution is that 
the value of n is not readily determined. Since A(5) = 0 at the 
boundary of the wake, say 5 = ±6, an infinite set of values of n would 
arise as roots of the equation, 


E 

i=l 


(2n-2i+2)(2n-2i+4) 

(2iTr~ 


(2n+2) 2ri 


0 


It would be necessary to test each value to find the best fit for the 
observations - a monumental taski 
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A third possibility arises by a change in definition of the 
quantity, £i(x). A somewhat unsatisfying feature of the previous 
solutions is that the width of the obstacle is not an explicit feature 
of the scaling. Let ?i(x) be redefined as L, the half-width of the 
cross section presented to the wind by the obstacle. L will vary as 
the direction of the mean reference wind varies. We assume as a 
solution, 

- U‘ = Uo(-QiC^- + Qi) A(x) . (32) 

where x = x/h, 5 = y/L, ana h ~ ^2 is the height of the structure. 

The quantity (-QiC^ + Qi) represents a parabola which requires U' 0 
at ^ = ±1 , i .e. , y = ±L. 

When (32) is substituted into (24), it is found that 

A(x) = , B = constant, 

and 

2ZkhU* 

" i7~ Uo(Z)(i-5^) 


It can be seen that 

U' e'^^^ 
" Uo(Z) 


QiBid-?^) . 


(33) 


Finally, a solution which actually differs very little from (32) 
is given by 


- U' = UoQ 2 Cos(|d A(x) (34) 

In this case it is found that 
A(x) = 


1 
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where 

Tr2ZkhU* 



■^2 - 4L2 Uo(zy * 


Then, 

- e^2X , Q2B2COs(f5) . 

(35) 


and since presumably U' = when C = 0, Q 2 B 2 can be determined from 
a plot of the left-hand side of (35) against 5 . 

Testing of Theory 

Table 3 contains a listing of those data which were suitable for 
testing the theoretical development. Most of these data were measured 
when the instruments were 8 ft above the ground; however, on 6-6-76 
the instruments were 15 ft above ground during runs 2 and 9, 30 ft 
above ground during runs 3 and 10, and 60 ft above ground during runs 
4 and 11. 

The third and fourth columns of Table 3 give the velocity com- 
ponents relative to a cartesian coordinate system in which x is in the 
direction of the reference wind (Uq) and passes through the center 
of the building; x = 0 was taken to be at the building wall on the 
downwind side. Columns 5 and 6 give the location of the observation 
point in this coordinate system. In this coordinate system observa- 
tions at points outside of y = ± L do not have a defined x position. 
Therefore, some data appearing in Tables 1 and 2 do not appear in 
Table 3. 

Experimental variability precluded a very good match of the 
observations with any of the theoretical relationships given previously. 
However, the most acceptable results were obtained using (35). 
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Table 3.’- Wake Data 


Date and 
Run 

Tower 

Relative Components 
U(m/sec) V(m/sec) 

Position 

X (ft) y (ft) 


1 

Uo 

= 3.0 

M 


4-17-74 

2 

2.4 

-0.2 

131 

-42.5 

Run 8 

3 

2.2 

0.2 

123 

- 3 


4 

2.1 

-0.7 

122 

36 


1 

Uo 

= 2.5 


_ _ _ 

4-23-74 

2 

1.7 

-0.5 

190 

-52 

Run 1 

3 

1.0 

0.7 

92 

1 


4 

1.0 

0.3 

92 

36 


5 

1.5 

-0.8 

92 

70.5 


4-24-74 
Run 1 

1 

5 

Uo = 
1.4 

2.0 

-0.5 

140 

-13 

10-16-75 

1 

Uo = 

= 6.8 


_ ^ _ 

Run 1 

4 

-2.0 

-2. 5 

45 

-45.4 

10-16-75 

1 

Uo = 

= 5.1 

„ ^ _ 

...... 

Run 2 

2 

0.0 

-1.3 

87 

-54.7 


4 

0.8 

-2.7 

87 

-17.8 

10-16-75 

1 

Uo 

= 3.6 

mm mm ^ 

_ _ __ 

Run 3 

2 

3.0 

-1.1 

146 

20.1 


4 

4.1 

-0.6 

258 

58.5 

10-17-75 

1 

Uo 

= 3.4 

B » av 

_ _ _ 

Run 1 

4 

-0.5 

-1.0 

45 

-42.9 


5 

-0.8 

-2.7 

45 

- 6.6 


10-17-75 
Run 2 

1 

2 

4 

Uo 

1.7 

3.8 

= 5.2 

-1.4 
0.2 

88 

167 

10.6 

46.7 

10-17-75 

1 

Uo 

= 5.4 

— 

— 

Run 3 

2 

1.1 

-1.9 

163 

-63.6 


4 

4.4 

-1.8 

163 

-34.2 
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Table 3. - Continued 


Date and 
Run 

Tower 

Relative Components 
U(m/sec) V(m/sec) 

Position 
X (ft) y (ft) 

10-17-75 

1 

Uo 

= 6.0 

iB as 


Run 5 

4 

4.0 

-2.2 

202 

-49.1 

10-17-75 

1 

Uo 

= 5.3 



Run 6 

2 

0.2 

-2,2 

77 

-45.3 


4 

3.4 

-2.3 

155 

-57.1 

6-6-76 

1 

Uo 

= 2.9 



Run 1 

2 

2.8 

0.2 

239 

47.8 

6-6-76 

1 

Uo 

= 3.5 


___ 

Run 8 

2 

2.7 

-0.4 

194 

-37.1 

6-6-76 

1 

Uo 

= 3.2 



Run 2 

2 

2.6 

-0.4 

247 

58 

6-6-76 

1 

Uo 

= 3.7 



Run 3 

2 

4.0 

-1.0 

278 

79 

6-6-76 

1 

Uo 

= 3.2 

~ W 


Run 4 

2 

3.5 

-0.7 

251 

42 

6-6-76 

1 

Uo 

= 3.6 


___ 

Run 9 

2 

2.6 

0 

179 

13 

6-6-76 

1 

Uo 

= 3.2 

___ 

^ _ 

Run 10 

2 

3.0 

0.2 

208 

-52 

6-6-76 

1 

Uo 

= 3.8 



Run 11 

2 

3.9 

-0.1 

195 

-25 
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Figure 5 shows a plot of - U'e^^/Uo(Z) vs. 5 , where 

. _ k^hTT^Z 

* " 4L^lnZ/zo ’ 

and zq was estimated to be 0.1 ft. The points shown as solid dots 
in Figure 5 are the only ones which are considered valid for purposes 
of testing the theory. These are the only points that satisfy the 
assumption of the theory that U' should be small compared to Uq. The 
curve ? = 0.4 cos(ir^/2) has been drawn as a suggested fit to these 
points - it has not been verified by the least-squares procedure. 

The reason the other points have been ignored also can be 
illustrated by reference to Figure 6 which shows wind vectors, 
normalized by division of the magnitude by Uo» at the various observa- 
tion points in the x-y coordinate system. This figure is only 
roughly representative of the flow that might be expected with 
simultaneous observations at many points. Obviously it does not 
account for variations in the shape of the cavity and the structure 
of the wake flow with variations in the orientation of Uq relative 
to the building, as shown by Evans (op. cit.). However, certain 
features in Figure 6 are indicative. For example, the points shown 
by crosses must represent observations within or very close to the 
cavity. The theory on the other hand was developed for the region 
downwind from the cavity. The observations with solid dots in 
Figure 6 correspond to those shown by solid dots in Figure 5, which 
indeed are downwind of the cavity region. 

The points shown by open squares in the two figures are observa- 
tions which seem to be in a transition region between the cavity 
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and the remainder of the wake flow. These observations also fail 
to meet the criterion that U' be small compared to Uq. The theory 
appears to become valid for values of x greater than 100 to 120 ft. 
This corresponds in the present case to values of 7 of 2.5 to 3. 

It should be pointed out that we were unsuccessful in obtaining 
usable observations at distances greater than seven building heights. 
This was because when working at these distances and moving the 
instruments about, the natural variability of the wind makes it 
difficult to guess where the boundaries of the wake will be in the 


mean. 
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CONCLUSIONS 

Based upon scaling and dimensional analysis, the statement 
-U' = Uo 6 cos(tt 5/2) e"*^^ , 

was derived for the flow in the wake of a building at points downstream 
of the cavity. In this statement, p is a constant with a value of 
about 0.4 and 6 depends upon the orientation of the undisturbed 
mean flow relative to the building. The vertical dimension of the 
structure is accounted for in x and the lateral dimension is approxi- 
mated in ?. 

The foregoing relationship appears to be a reasonable fit to the 
observations; however, there is a large degree of experimental varia- 
bility and additional observations, especially at distances greater 
than seven building heights, are needed to provide complete verification. 

The results indicate that an aircraft passing through the wake 
during take-off or landing would experience not only a change in 
turbulence level, but a change in mean wind speed of a magnitude 
roughly equivalent to that of the eddy components. 

An interesting difference between the results of wind tunnel 
studies and those obtained here is with regard to the bubble or 
cavity. Depending upon the author, reattachment has been found in model 
studies to occur seven or more obstacle heights downstream. Evans' 

(op. cit.) results are the exception. For a building of the size and 
shape of the hangar involved in this study, his results indicate that 
reattachment should have occurred about two and a half building heights 
downstream. Our results are in good agreement with this. 
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